DEVELOPMENT OF A CONCEPTUAL MODEL AND PARAMETERIZATION OF A MATHEMATICAL MODEL FOR SEDIMENT PHOSPHORUS DIAGENESIS, SORPTION AND RELEASE FROM ONONDAGA LAKE SYRACUSE, NEW YORK

Eutrophication is a persistent problem in many fresh water lakes. Delay in lake recovery following reductions in external loading of phosphorus, the limiting nutrient in fresh water ecosystems, is often observed. Models have been created to assist with lake remediation efforts, however, the application of management tools to sediment diagenesis is often neglected due to conceptual and mathematical complexity. SED2K (Chapra et al. 2012) is proposed as a "middle way", offering engineering rigor while being accessible to users. An objective of this research is to further support the development and application SED2K for sediment phosphorus diagenesis and release to the water column of Onondaga Lake. Application of SED2K has been made to eutrophic Lake Alice in Minnesota. The more homogenous sediment characteristics of Lake Alice, compared with the industrially polluted sediment layers of Onondaga Lake, allowed for an invariant rate coefficient to be applied to describe first order decay kinetics of phosphorus. When a similar approach was attempted on Onondaga Lake an invariant rate coefficient failed to simulate the sediment phosphorus profile. Therefore, labile P was accounted for by progressive preservation after burial and a rate coefficient which gradual decreased with depth was applied. In this study, profile sediment samples were chemically extracted into five operationally-defined fractions: CaCO3-P, Fe/Al-P, Biogenic-P, Ca Mineral-P and Residual-P. Chemical fractionation data, from this study, showed that preservation is not the only mechanism by which phosphorus may be maintained in a non-reactive state in the profile. Sorption has been shown to contribute substantially to P burial within the profile. A new kinetic approach involving partitioning of P into process based fractions is applied here. Results from this approach indicate that labile P (Ca Mineral and Organic P) is contributing to internal P loading to Onondaga Lake, through diagenesis and diffusion to the water column, while the sorbed P fraction (Fe/Al-P and CaCO3-P) is remaining consistent. Sediment profile concentrations of labile and total phosphorus at time of deposition were also modeled and compared with current labile and total phosphorus, to quantify the extent to which remaining phosphorus which will continue to contribute to internal P loading and influence the trophic status of Onondaga Lake. Results presented here also allowed for estimation of the depth of the active sediment layer and the attendant response time as well as the sediment burden of labile P and associated efflux.

Hydrogeochemical controls on uranium in aquifers of the Jacobsville Sandstone

The occurrence of elevated uranium (U) in sandstone aquifers was investigated in the Upper Peninsula of Michigan, focusing on aquifers of the Jacobsville Sandstone. The hydrogeochemical controls on groundwater U concentrations were characterized using a combination of water sampling and spectral gamma-ray logging of sandstone cliffs and residential water wells. 235U/238U isotope ratios were consistent with naturally occurring U. Approximately 25% of the 270 wells tested had U concentrations above the U.S. Environmental Protection Agency Maximum Contaminant Level (MCL) of 30 μg/L, with elevated U generally occurring in localized clusters. Water wells were logged to determine whether groundwater U anomalies could be explained by the heterogeneous distribution of U in the sandstone. Not all wells with relative U enrichment in the sandstone produced water with U above the MCL, indicating that the effect of U enrichment in the sandstone may be modified by other hydrogeochemical factors. Well water had high redox, indicating U is in its highly soluble (VI) valence. Equilibrium modeling indicated that aqueous U is complexed with carbonates. In general, wells with elevated U concentrations had low 235U/238U activity ratios. However, in some areas U concentrations and 235U/238U activity ratios were simultaneously high, possibly indicating differences in rock-water interactions. Limited groundwater age dating suggested that residence time may also help explain variations in well water U concentrations. Low levels of U enrichment (4 to 30 ppm) in the Jacobsville sandstone may make wells in its oxidized aquifers at risk for U concentrations above the MCL. On average, U concentrations were highest in heavy mineral and clay layers and rip up conglomerates. Uranium concentrations above 4 ppm also occurred in siltstones, sandstones and conglomerates. Uranium enrichment was likely controlled by deposition processes, sorption to clays, and groundwater flow, which was controlled by permeability variations in the sandstone. Low levels of U enrichment were found at deltaic, lacustrine and alluvial fan deposits and were not isolated to specific depositional environments.

Simulating dioxane transport in a heterogeneous glacial aquifer system (Washtenaw County, Michigan) using publicly available models and data

The primary challenge in groundwater and contaminant transport modeling is obtaining the data needed for constructing, calibrating and testing the models. Large amounts of data are necessary for describing the hydrostratigraphy in areas with complex geology. Increasingly states are making spatial data available that can be used for input to groundwater flow models. The appropriateness of this data for large-scale flow systems has not been tested. This study focuses on modeling a plume of 1,4-dioxane in a heterogeneous aquifer system in Scio Township, Washtenaw County, Michigan. The analysis consisted of: (1) characterization of hydrogeology of the area and construction of a conceptual model based on publicly available spatial data, (2) development and calibration of a regional flow model for the site, (3) conversion of the regional model to a more highly resolved local model, (4) simulation of the dioxane plume, and (5) evaluation of the model's ability to simulate field data and estimation of the possible dioxane sources and subsequent migration until maximum concentrations are at or below the Michigan Department of Environmental Quality's residential cleanup standard for groundwater (85 ppb). MODFLOW-2000 and MT3D programs were utilized to simulate the groundwater flow and the development and movement of the 1, 4-dioxane plume, respectively. MODFLOW simulates transient groundwater flow in a quasi-3-dimensional sense, subject to a variety of boundary conditions that can simulate recharge, pumping, and surface-/groundwater interactions. MT3D simulates solute advection with groundwater flow (using the flow solution from MODFLOW), dispersion, source/sink mixing, and chemical reaction of contaminants. This modeling approach was successful at simulating the groundwater flows by calibrating recharge and hydraulic conductivities. The plume transport was adequately simulated using literature dispersivity and sorption coefficients, although the plume geometries were not well constrained.